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O ' ABSTRACT 

o ' 



We consider the effects of inverse Compton scattering of synchrotron photons from 
relativistic electrons in GRB afterglows. We compute the spectrum of the inverse Comp- 
ton emission and find that it can dominate the total cooling rate of the afterglow for 
several months or even years after the initial explosion. We demonstrate that the pres- 
ence of strong inverse Compton cooling can be deduced from the effect it has on the 
time-evolution of the cooling break in the synchrotron spectral component, and therefore 
on the optical and X-ray afterglow lightcurves. We then show how the physical inter- 
im [ pretation of the observed characteristics of the synchrotron spectrum must be modified 
to take into consideration this extra source of cooling, and give a revised prescription 
for computing physical parameters characterizing the expanding shock wave from the 
observed quantities. We find that for a given set of observables (synchrotron break 
frequencies and fluxes) there is either no consistent physical interpretation or two of 
them. Finally we discuss the prospects of directly detecting the inverse Compton emis- 
sion with Chandra. We argue that such a detection is possible for GRBs exploding in 
£j | a reasonably dense (n > 1cm -3 ) medium. 
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1. Introduction 

It is widely accepted that the emission from GRB afterglows is produced in relativistic shocks, 
as the ejecta from an underlying explosion expands into the surrounding medium (e.g. Piran 1999 
and references therein) . In the standard picture, the ambient electrons are accelerated in the shock 
front to highly relativistic energies, with their Lorenz factors described by a power law distribution 
above some minimum value. Besides particle acceleration, the shock is also responsible for the 
creation of reasonably strong magnetic fields. Under these conditions synchrotron radiation is 
produced, whose spectrum and light curves was described by Sari, Piran & Narayan (1998). 

Though GRBs and their afterglow are optically thin to electron scattering for years after 
the initial explosion, some synchrotron photons will Compton scatter on the shock-accelerated 
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electrons, producing an additional inverse Compton component at higher energies (Panaitescu & 
Meszaros 1998; Wei & Lu 1998; Totani 1998; Chiang & Dernier 1999; Panaitescu k Kumar 2000; 
similar calculations in other astrophysical contexts were done by Blandford & McKee 1977 both 
for a constant density and a wind- like ambient medium). Since only a negligible fraction of the 
synchrotron photons will be scattered, this will have no direct effect on the shape of the synchrotron 
spectrum. However, the ratio of inverse Compton to synchrotron luminosities is equal to the square 
root of the ratio of the electron and magnetic energy densities behind the shock front (Sari, Narayan 
& Piran 1996). Since this number can be significantly above unity, this implies that the electron 
cooling rate via inverse Compton scattering of synchrotron photons must be considered in order to 
create a realistic physical description of the GRB afterglow emission. 

In their paper, Panaitescu & Kumar ( 2000) give a detailed treatment of the inverse Compton 
emission from GRB afterglows. Here we extend their analysis, concentrating on the implications 
of this cooling process for the observable properties of these objects. We start by calculating the 
spectrum of the inverse Compton emission in §2. The gross power-law features of this spectral 
component are the same as those in Panaitescu & Kumar ( 2000). However, here we show that the 
broken power-law is not an adequate approximation above the peak frequency, where the presence 
of the logarithmic terms may introduce large corrections to the predicted flux levels. These extra 
terms also serve to smooth out the resulting spectrum and high-energy lightcurves, eliminating 
sharp power-law breaks in the theoretical predictions. 

In §3 we discuss in detail the effects of inverse Compton cooling on the afterglow evolution. We 
show that instead of the two standard stages of evolution, corresponding to fast and slow cooling, 
we now have three stages. This occurs because the slow cooling regime should now be divided 
into two, the early stage during which inverse Compton scattering dominates the total cooling, and 
the late stage when it is unimportant compared to energy loss via synchrotron. In §4 we focus on 
how the physical interpretation of the observed spectral break frequencies and the time evolution 
of the afterglow emission must be modified to take into account the presence of extra cooling via 
inverse Compton scattering. In particular, we derive a useful physical limit on a combination of 
the observable parameters, measurements of which are generally the goal of afterglow observations. 
This result is then discussed in the context of GRB 970508. 

Finally, in §5 we discuss the possibility of observing the IC emission directly. We show that the 
prospects of detecting this spectral component are very good for GRBs which explode in a medium 
with density of order 1 cm~ 3 . Our conclusions are summarized in §6. 

2. Spectrum 

We make the standard assumption of a uniform magnetic field and a power-law injection of 
electrons, with the energy distribution given by N(j) oc 7~ p , behind the expanding shock front. 
Then the general shape of the afterglow spectrum is determined by the relationship between the 
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electron cooling time and the system dynamical time. Below we consider two limiting cases: when 
the majority of the injected electrons can cool on the dynamical time of the system (fast cooling) or 
when the cooling affects only the electrons in the high energy tail of the distribution (slow cooling) . 



In this case the electron energy distribution around the minimum injection energy, j m , is not 
affected by the cooling. However, electrons with Lorenz factors above some critical value, 7 C , radiate 
significant fraction of their energy on the dynamical time. Since the ratio of the cooling time of the 
electron with Lorenz factor 7 e to the dynamical time is (7 c /7e)> the electron distribution above 7 C 
become steeper by one power. Therefore, the resulting electron energy distribution is : 



As described by Sari et al. (1998), the synchrotron emission from such distribution of electrons 
can be approximated by a broken power-law spectrum with three characteristic break frequencies. 
One is the self-absorption frequency, v a , below which the system becomes optically thick. The 
other two are the peak frequencies of the emission from the electrons with the characteristic Lorenz 
factors 7 m and 7 C , denoted v m and v c , respectively. Here we assume that v m and v c are greater 
than v a . At frequencies v < v a , self-absorption is important and the specific flux is proportional 
to u 2 , while between v a and v m the synchrotron emission grows more slowly with / oc v 1 / 3 . The 
specific flux peaks at v m and then decreases as u^^ 1 )/ 2 for frequencies v m < v < v c \ above v c it 
decreases as v~ v l 2 . 

Using this synchrotron spectrum as the source of seed photons, we computed the resulting 
inverse Compton emission, by integrating over the differential scattering cross section and over the 
electron energy distribution given by Eq. 2-1 (this calculation is discussed in detail in appendix A). 
We find that, like the original synchrotron spectrum, the upscattered component consists of four 
segments, with the breaks at three characteristic frequencies: v]p ~ 27 m z/ a , ~ 27 m z/ m , and 
v]F ~ 2"fc u c- At frequencies v < u^, the main contribution to the IC spectral component comes 
from synchrotron photons scattered by the electrons with Lorenz factor j m . In this region, the 
spectrum is well approximated by two power-law segments. The specific flux initially increases 
linearly with frequency up to v[ ; this spectral slope is determined by the shape of the scattering 
cross section (see appendix A). It then continues to rise as v x l z up to its peak at v = v 1 ^ . 

In the frequency range vj~ < u < v~ , electrons with a range of Lorenz factors between 7 m and 
7 C contribute equally to the emission at each frequency. Because of this, the flux drops off roughly 
as z/^p -1 )/ 2 , but with two major differences from the original synchrotron spectrum. Firstly, the 
region with this spectral slope extends over double the frequency range (in logarithmic units) of 
the corresponding synchrotron component, i.e. ' jv^F = (y m jv c ) 2 . Secondly, contributions from 



2.1. Slow Cooling 




7m < 7 < 7c 

7c < 7- 



(2-1) 
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different electron energies generate a logarithmic term, which peaks at v = \J v^vl c , increasing 
the flux there by a factor of order In (v c /u m ) above that of the the underlying power-law (see Eq. 
[AS]). 

Similarly, at frequencies above v]P ', the emission is composed of equal contributions from 
electrons with a range of Lorenz factors 7 e > 7 C . Here specific flux drops off as v~' p l 2 , but also with 
an addition of a logarithmically growing term of order In iy jv], ) (see Eq. [A5]). 

The inverse Compton spectral component, approximated as a broken power-law and normalized 
using Eq. (2-3), below is shown in Fig. la together with the more "exact" spectrum obtained by 
computing the integrals in Eq. (A2). To emphasize the importance of IC emission as a cooling 
mechanism, we plotted vf v vs. v, which gives a measure of energy emitted per logarithmic frequency 
interval. The comparison of the two curves shows that the power-law approximation does a very 
good job at frequencies below , but becomes somewhat inadequate at higher frequencies. In 
fact, our calculation shows that at frequencies v > ', the spectrum of the IC emission has a 
continuously varying slope, without identifiable spectral breaks. Moreover, one must keep in mind 
the the observed IC spectral component is expected to be even smoother, since in our calculations 
we assumed a broken power-law distribution of the synchrotron seed photons and used a simplified 
description of the scattering cross section (see appendix A). 

Another feature of the detailed spectrum, which cannot be deduced from the power-law ap- 
proximation, is that for small values of p the IC energy output actually peaks well above is[ c (which 
is the peak emission frequency for the approximate power-law spectrum). It is easy to show that 
this feature is due to the presence of the logarithmic correction term. It persists for 2 < p < \/6, 
which includes the range of values deduced from observations (e.g. Wijers & Galama 1999; Sari, 
Piran &, Halpern 1999; Freedman h Waxman 2000). 



2.2. Fast Cooling 

In this regime all of the injected electrons are able to cool on the dynamical time scale of the 
afterglow. Therefore, there is a population of electrons with Lorenz factors below the injection 
minimum 7 m , and the resulting electron distribution takes the form 

AT( 7 ) oc { 7 " 2 ' , * 7c<7<7m5 (2-2) 
UJ \ 7 > if lm < 7- K ' 

The synchrotron spectrum has the same three characteristic break frequencies, though now 
v c < v m , and consists of four power-law segments. Below v a the emission is optically thick and the 
specific flux increases as u 2 , and in the range u a < v < v c we have f v oc v 1 ^. The flux peaks at v c 
and declines as / oc u~ 1 / 2 for v c < v < v m and as / oc v~ p l 2 above u m . 



A spectrum of the inverse Compton emission in this regime is plotted in Fig. lb. Its general 
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features are very similar to those of the slow cooling case described above. The specific flux is linear 
in frequency below v]p and continues as / oc v 1 ^ up to its peak at u[ c , in the range v]F < v < z/^p 
it declines as / oc i/ -1 / 2 and drops off with / oc v~ p l 2 above ■ Since most of the contribution to 
the IC emission below v~ comes from the lowest energy electrons, a broken power-law is a fairly 
good approximation to the spectrum in that region. At higher frequencies electrons with a range 
of Lorenz factors contribute to the emission, creating additional logarithmic terms (see Eq. [A8]), 
which smooth out the breaks in the spectral slope. 



2.3. Peak Flux and Luminosity Ratios 

The strength of the inverse Compton emission relative to synchrotron can be estimated by 
considering the ratio of specific fluxes measured at the peak of the respective spectral components 
( peak of f v rather than vf v ). We denote these peak fluxes as f ma x an d fmax f° r the synchrotron 
and inverse Compton components, respectively. As described in §2.1 and §2.2 above, the former 
peaks at v m and the latter at in the slow cooling regime; and at v c and v£ c , respectively, in the 
fast cooling regime. In both cases the number of electrons that contribute to the radiation around 
the maximum is proportional to the total number of electrons in the system, N. Therefore the 
ratio of the two fluxes is simply 

~ T^2 ~ I^tuR = 2 x 10- 7 ni R 18 , (2-3) 

J max 47Trt o 

where the last two terms assume that the electrons have been collected from an ambient medium 
with particle density n over a distance R. For a range of values of n and R, which are relevant to 
the observable GRB afterglows, the flux ratio above is always much less than unity. 

To assess the relative importance of contributions from the inverse Compton and synchrotron 
emission to the total cooling rate of the electrons, we can compute the ratio of the total energies 
emitted via these two mechanisms. Instead of integrating over the entire spectrum, we evaluate 
vj v at the peaks of the two spectral components and simply compute a ratio of the two values. 
Note that taking L syn ~ v c f v (v c ) and Lie ~ v c C fv{ v c C ) is not m general a good assumption, since 
direct spectral integration gives luminosities a factor >10 greater than this estimate. However, 
since the shapes of the synchrotron and IC spectral components are similar, the correction factors 
nearly cancel out. 

In the slow cooling regime, the energy emitted via synchrotron peaks at v c . The inverse 
Compton emission can reach its maximum energy output at v 3> v]P for small values of p (see 
§2.1), however, since vj v is nearly flat in this region of the spectrum, for the purposes of this 
estimate we use v\ fviyi )- Then the luminosity ratio during the slow cooling regime is 

^-Wf^) 1 " (2-4) 
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In the fast cooling regime the synchrotron and the IC energy emission peak at u m and , 
respectively, and the luminosity ratio is given by 

~ |cr T ni?7c7 m . (2-5) 

-'-'syn " 



3. Afterglow Evolution with Strong Compton Cooling 

Throughout most of the early afterglow evolution, inverse Compton emission typically domi- 
nates the total cooling rate, and therefore it determines the value of break energy in the electron 
energy distribution. This has a significant effect on the time evolution of both the synchrotron and 
the inverse Compton component. 



3.1. Luminosity ratio 

The ratio of the inverse Compton to synchrotron luminosity can be computed in a more 
general way (Sari, Narayan & Piran 1996), that does not deal with the details of the spectrum, but 
depends only on the underlying physical properties of the expanding shock wave . We generalize the 
derivation given by Sari et al. (1996) to describe both fast and slow cooling regimes by introducing 
a parameter rj, equal to the fraction of the electron energy that was radiated away (via both 
synchrotron and IC emission). Then the ratio of luminosities, in the limit of single scattering, is 
given by 

= L IC = Urad = Usyn = T]U e /{l+x) ^ Tje e 

X ~L syn ~ U B ~ U B U B ~e B (l + xY { ' } 

where U syn , U B and U e are the energy density of synchrotron radiation, magnetic field and rela- 
tivistic electrons, respectively. Note that in general U syn = r](3U e /(l + x), where (3 is the velocity 
of material behind the shock front (in the frame of the shock); however, for a relativistic shock (3 is 
of order unity. The parameters e e and e b are defined as fractions of the total explosion energy that 
go into accelerating electrons and amplifying magnetic fields behind the shock front, respectively. 

Solving Eq. (3-1) for x we obtain 



-i+ ji+^a 

x = V- B -. (3-2) 



This solution has two interesting limits: 



if 2££<1, 



1 (5) 



Hereafter, unless specified otherwise, we always use these approximate solutions for x. 
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Clearly, if ??e e /e j B <C 1, the inverse Compton cooling rate is unimportant compared to syn- 
chrotron, and if r]e e /eB 3> 1, it dominates the total emission. Using these expressions, the relative 
importance of IC emission can be evaluated using only the fundamental properties of the expand- 
ing shock, and the ability of the electrons to cool. We can reduce this result to that in Eqs. (2-4, 
2-5) above by noting that (i) for slow cooling rj = ('Jc/lm) 2 ~ p and for fast cooling ij = 1; (ii) 
7m ~ e e U/m e c 2 n, where U is the total energy density behind the shock; (hi) j c ~ m e c 2 / aTcesUt 
if rje e /eB <C 1, and 7 C ~ m e c 2 /xaTcesUt if rje e /eB S> 1. 

3.2. Effects on Synchrotron Spectrum and Lightcurves 

Inverse Compton scattering can, in principle, affect the observed spectrum in three different 
ways (see also discussion in Sari, Narayan and Piran 1996). In the first place, it reduces the num- 
ber of seed photons, changing the overall normalization of the synchrotron spectral component. 
This is important only if the Thompson optical depth is larger than unity, and is therefore neg- 
ligible in GRB afterglows. Secondly, IC scattering produces an additional emission component at 
higher frequencies, which we described it §2. Finally, when the inverse Compton emission domi- 
nates the overall cooling of the electrons, it reduces the energy available for synchrotron radiation. 
Consequently, the cooling break energy of the electron distribution, j c , is reduced from its stan- 
dard (synchrotron only) value by a factor (1 + x), i.e. in all our calculations we must replace 
7c — ► 7c/(l + x )- It is the results of this latter effect that we discuss below. 

From Eq. 3-3 it is clear that if e e < e#, the IC power is never larger than the synchrotron 
power, since by definition rj < 1. However, current estimates based on afterglow observations, show 
that es < e e ; in fact in several objects the ratio e e /eB seems to be close to a hundred (Wijers &: 
Galama 1999; Granot et al. 1999). This implies that the magnetic field is relatively weak, and IC 
emission is important as long as rf > 0.01. 

Assuming that e e > e#, the evolution of a GRB afterglow must now be divided into three 
regimes. In the first (fast cooling) stage, most of the electrons energy is lost, so that r/ = 1 and 
IC emission dominates over synchrotron by a time-constant factor y / e e /e_e. Later, during the slow 
cooling period, r\ decreases and so is the importance of IC scattering, though it still dominates the 
total cooling. Finally, r/ becomes so small that only synchrotron emission is important; this is the 
traditional slow cooling regime. Below we discuss each of these stages in detail. Note that all our 
numerical results are computed in terms of es-2 = es/0.01 and e ej .5 = e e /0.5. 

Since temporal evolution of the afterglow emission depends on the external density profile, the 
results are different in the case of the uniform ambient medium and the wind scenario, in which the 
density decreases as p oc r~ 2 . In this section we will present the results for the constant ambient 
density, and discuss the wind-like density profile in appendix B. 
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3.2.1. Fast cooling 

During the fast cooling stage, IC scattering increases the radiation losses of each electron by 
a factor 1 + x ~ x ~ y^T/^B- Therefore, the cooling frequency v c oc 7 2 is reduced from its 
synchrotron-only value by a factor x 2 ~ e e /e£. Moreover, since low energy synchrotron emission 
is dominated by the electrons near the cooling break, the synchrotron self-absorption frequency 
increases by a factor x ~ y/ee/cB- Because both e e and es remain constant in time, these changes 
do not affect the time evolution of the synchrotron spectrum. 

Since IC scattering increases the total cooling rate, it prolongs the duration of the fast cooling 
regime. For a constant ambient density the transition to slow cooling is delayed by a factor e e /e£, 
if e e > €b- The correct transition time then becomes 

ti c = 170 (1 + z)e 3 e e B E 52 n 1 days = 6.3 (1 + z)e 3 . 5 e B -2#52ni hours, (3-4) 

where E = E52 x 10 52 ergs _1 is the total energy of the explosion, and n = ri\ x 1 cm~ 3 is the 
ambient electron density. By comparison, if e e /eB < 1, IC scattering is never an efficient cooling 
mechanism and we get the usual (synchrotron only) expression: 

t = 170 (1 + z)e 2 e e 2 B E 52 n 1 days. (3-5) 



3.2.2. Slow IC-Dominated Cooling 



During the slow cooling stage only the cooling frequency is affected by IC cooling. The pa- 
rameter 77 is no longer equal to unity, so during this intermediate regime v c is reduced by a factor 



(l+x) 



x 



r]e e /eB, as long as x > 1. It is important to note that though 77 decreases with time, 



its instantaneous value can be estimated directly from the observed synchrotron spectrum alone, 

p — 2 

since in the slow cooling regime 77 = (7 c /7m) 2 ~ p = {yd v m) ~ ■ Since the value of p is estimated 
to be in the range of 2.2 — 2.4 the decrease of 77 is very slow and inverse Compton scattering remains 
a dominant cooling process for a long time. 

Without inverse Compton cooling, v c oc t~ x l 2 and v m oc t~ 3 / 2 so that their ratio is v c jv m = t/to. 
Taking IC emission into account we have v c /is m ~ (t/to)x~ 2 , and x evolves with time according to 
the following equation, 



which gives 



x = 



x ~ 





1 

(4-p) 



(P-2) 
2(4-p) 




(P-2) 
2(4-p) 



(3-6) 



(3-7) 



One should keep in mind however, that in deriving the equation above we used the relativistic 
time scalings for the frequency ratio u c /u m . Once the fireball is no longer relativistic, our basic 
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expression for x (Eq. 3-1) is no longer valid. In the non-relativistic case x decreases as i~ 3 / 10 for 
p = 2, and even more rapidly for larger p, considerably faster than implied in Eq. (3-7). Therefore, 
IC cooling ceases to be dominant shortly after the relativistic stage is over. 



3.2.3. Slow Synchrotron- Dominated Cooling 

The third stage of afterglow evolution, when x < 1 and the IC cooling rate is weaker than the 
synchrotron cooling rate and therefore, has no effect on the synchrotron spectrum, begins at time 

4—p 4—p 

,IC _ AC ( e e \ <^ _ ( Ce,.5 \ ~ 2 ,3 



\7^) =3yCarS Vef^J e t^B,-2E 52 ni, (3-8) 

with the numerical coefficient computed for p = 2.3. Note that this expression is very sensitive to 
the value of p, and goes to infinity when p — > 2. 

Equation 3-8 implies that for e e /e# > 10, the IC scattering dominates total cooling over the 
whole relativistic stage of afterglow evolution, and therefore, over the time period containing most 
current observations. 



3.3. Inverse Compton Lightcurves 

Equations describing light curves of inverse Compton emission in the regimes of fast cooling 
and slow cooling with negligible IC (as defined in §3.2.1 and §3.2.3) were given in Panaitescu & 
Kumar ( 2000), so we omit them here. One must keep in mind however, that these equations were 
derived using the broken power-law approximation to the IC spectral component, and so must be 
used with caution. During the intermediate stage (see §3.2.2), when v c > v, m but the total cooling 
rate is dominated by IC scattering, the break frequencies of the IC spectral component and its peak 
flux evolve with time in the following way 

v{ c = 2 7 2 m » a xt- 3 /\ (3-9) 

= 2 7 ^ m «r 9 / 4 , 

vf = 2 7 > c oc(l + x)- 4 t- 1 /4 ocr 1 / 4+ ^, 
fi c ~ cr T Rnf m oc t x '\ 

For a range of values of p consistent with constraints from observations, the break frequency v]F 
remains nearly constant or increases with time: v]P oc t~ - 03 for p = 2.2 and v]F oc t ' 25 for p = 2.4. 
It is therefore unlikely that it can ever be observed directly with X-ray instruments (see Fig. 1). 

Using the expressions above and the detailed spectrum described in the appendix, one can 
obtain the lightcurve of the IC emission at any frequency. As we have emphasized earlier, taking 
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V < V, 
,,IC 



IC 



only the broken power-law approximation to the IC spectrum can miss considerably, and can be 
used only to get a very rough estimate, given below 

( t 9 / 4 , if 
t\ if 
is, if 

9p-10 , p-2 

t 8 +4 -p, if 



< V < V, 



IC 



U 1 ^ < V < v 1 ^ 
v{ c < V. 



(3-10) 



The equation above differs from that given by Panaitescu and Kumar above the frequency 
v]P . Their expression is valid only it late times, when IC cooling is not important, while Eq. (3- 
10) is valid during the intermediate stage of afterglow evolution (slow cooling dominated by IC). 
As we stressed earlier, for typical GRB parameters, this intermediate stage lasts throughout the 
relativistic evolution of the fireball. 



4. Inferring Afterglow Parameters from the Synchrotron Spectral Component 

Even if the IC emission is not directly observed, its presence affects how we interpret the 
properties of the observable synchrotron spectral component. The general shape of the synchrotron 
emission is not affected by the IC scattering. However, when this extra source of cooling is taken 
into account, it lowers the predicted value of the cooling frequency, v c by a factor (1 + x) 2 . In 
addition to this, since during the fast cooling stage the self absorption frequency, is a , depends on 
the the cooling frequency, in this regime its value increases by a factor (1 + x). The resulting 
modified expressions for the break frequencies and the peak flux of the synchrotron spectrum are 
listed below: 

vt° w = 3.6GHz(l + .)- 1 e -i. 5 4 /5 _ 2 i? 5 1 2 /5 nf , (4-1) 
v?" = 0.15 GHz (l + zr^e^EU'^t^l + x), (4-2) 
v m = bxW^Rzil + zy/hf^el^El^ 2 , (4-3) 
v c = 2.7 x 10 15 Hz (1 + z)-^e- B 3 LlE^ 2 n^t-^ 2 (l + x)" 2 , (4-4) 
fmax = 2.6mJy(l + z)e^_ 2 E 52 n 1 / D L28 , (4-5) 

where -Dl,28 is the luminosity distance in units of 10 28 cm. The coefficients in these equations were 
taken from the synchrotron spectrum calculation of Granot et al. (1999) and are slightly dependent 
on p. The numerical values here are quoted for p = 2.2. 

Combining Eqs. (3-1), (4-l)-(4-5), we can now derive an interesting constraint on the observ- 
able afterglow parameters: 

c . om+ ,r^#0 (c^) ¥ (j^f (W^) 1 (^)~ ' = jT^p- («) 
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The expression on the right hand side has a maximum value of 1/4. Thus, Eq. (4-6) above gives us 
a theoretical constraint on this combination of the observed break frequencies and peak flux for any 
afterglow spectrum, completely independent of its underlying physical parameters! Note that rj is 
also an observable quantity, since it is given by rj = min[(f c /z/ m )( 2-p )/ 2 , 1] The only assumption 
made in deriving this constraint is that we are seeing the synchrotron emission from a power-law 
distribution of electrons accelerated by a relativistic shock wave. 

For fast cooling, we can write a similar combination. It is simpler, however, to use same 
combination C as above, in which we substitute the "adjusted" frequency given by 

/ \ V 2 

slow fast I v c 



m 



<° w = < st — , (4-7) 



instead of the observed self absorption frequency u[ a st . 

To illustrate how this constraint can be used, let's consider GRB 970508, which has probably 
the best studied afterglow to date. Galama et al. (1998) described the observed broad band spec- 
trum of that burst and have shown it to be in good agreement with the synchrotron spectrum of 
Sari et al. (1998). Their fit shows that at t = 12.1 days u c ^ 1.6 x 10 14 Hz v m ^ 8.6 x 10 10 Hz, 
v a = 3.1GHz and f max = 1.7 mJy and p = 2.2. Using the break frequencies of GRB 970508 as 
given by Galama et al. (1998), we see that Eq. (4-6) gives C = 4 > 1/4. Though this result 
seems to suggest that the observed spectrum could not be produced by synchrotron emission, a 
more likely explanation is that the values of the break frequencies ax© not axxurcitG. Both Vra 

and 

v c are rather weakly constrained by the observations, and a factor of a few could easily resolve this 
disagreement. 

We recommend the following procedure to determine the afterglow parameters from a snapshot 
of the spectrum. First, one calculates the combinations of observable as given in Eq. 4-6 (with the 
correction of Eq. 4-7, if v c < u a ). If this combination of observable C is above 1/4, then there is 
no consistent solution. If C < 1/4 there are two solutions which can be found in the following way. 
First we solve Eq. (4-6) to give the two possible values of x: 



^ '-^ (4-8) 

and 

1-2C + V1-4C _ 1 , An . 

x 2 = ^ ^ > 1. (4-9) 

The two solutions for the physical parameters describing the afterglow are then given by: 

£ - = °' 23 (ds) " f (i^a) "* Go*e) " (^) § + 2) ~ 2 ^ s(1 + x)i (4 - 10) 

" ^ 23 ( es) § (iife) ^ (l^Hi) ' (^) & + 2)1D - (1 + ^ (4 -"> 
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* - « (AT 1 (lifer* W (&)*<u + + *>-» 

- - »» 72 (& * (loSk) 8 (life) H^)" 1 + 2)5 °- (1 + x,i (4 - 13) 

Since 1 + x\ = 1 one of these solutions degenerates to that given in Wijers & Galama (1999) 
with the corrected coefficients of Granot et al. (1999). However, in the second case we have 
X2 > 1 so that the resulting values of £52, c e , £b and n\ are significantly different. Substituting 
1 + X2 ~ X2 ~ 1/C in these equations we obtain 




This second solution where IC dominates, was neglected so far. It has a somewhat higher energy, 
higher density, higher e e and lower e#, than the low-x solution. 

How can we determine which of the two solutions is correct at a given moment? This cannot 
be done by studying a single snapshot of the afterglow spectrum. However, the temporal evolution 
above the cooling frequency is a little different in the two cases. The decay rate is slower in 
the second, IC-dominated, solution. If the decay rate below and above the cooling frequency is 
measured to be t~ ai and t~ a2 , respectively, then (for constant density of ambient material) the 
standard model, which ignores IC emission, predicts 02 — ol\ = 1/4. The presence of strong IC 
cooling will make the difference smaller, giving 012 — ol\ = 1/4 — (p — 2)/2(4 — p) = 0.2. For a 
wind- like ambient density profile (see appendix B), the standard model predicts «2 — ol\ = —1/4, 
which changes to ct2 — a\ = —1/4 — (p — 2)/(4 — p) = —0.36 when the IC emission dominates the 
cooling. This difference is not large (except perhaps in the wind model), and therefore, accurate 
measurement with long temporal baselines are needed. 

It was recently argued by Freedman & Waxman (2000) that the total energy contained in 
the electrons behind the shock e e E can be estimated using a single observation above the cooling 
frequency. They argued that this is insensitive to the value of the magnetic field. Their analysis, 
which ignores IC, could be explained as follows. Since the electron energy distribution is flat, i.e. 
p ~ 2, each decade in the electron distribution contains an amount of energy of the order of the 
total energy in all the electrons. Above the cooling frequency all of this energy is radiated and 
therefore, the observed vF v gives an approximate measure of the total electron energy. However, 
we have shown that if the ratio (e e /es) is large, most of this energy is radiated as IC emission 
rather than synchrotron. Therefore, during the fast cooling stage the energy estimate suggested by 
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Freedman & Waxman (2000) falls short of e e E by a factor \/e e /e£, while during the intermediate 
stage it is short by a factor ^/rfe^JeB, which is roughly proportional to e^ ' 45 ; and both of these 
correction factors depend on the magnetic field. Thus, Freedman & Waxman (2000) method can 
provide only a lower limit on e e E. 



5. Direct Detection of IC Emission 



Clearly, the IC emission component can be detected directly only at frequencies where it 
dominates over the synchrotron emission. As can be seen from Fig. 1, for typical parameters this 
occurs in the X-ray and 7-ray bands. Here we limit our discussion to the intermediate regime 
(defined in §3.2.2) of afterglow evolution, during which the crossing point between the synchrotron 
and IC spectral components, v IC , generally lies above the synchrotron cooling frequency u c and 
between IC break frequencies v]p and v~ . Using the power-law approximation for both spectral 
components it is then straight-forward to compute v IC ', both at early times, when v lc < ', and 



at later times, when v IC > v£ ■ We obtain the following expressions: 



.IC 



.IC 



,ic 



7c 
1m 



(27m 7c) 
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3 (3pf-8p-12) 
(X f2 (2+3p)(4-p) 



„ 3pf-23p+36 
^ C ^(2 7m7c )2- P OC t . 



if u IC < vl9- 



if u IC > vl9 . 



(5-1) 



7m 



2(p-2) 



7m 



It is easy to show that the transition between these two regimes occurs when 
Note that it is not necessary to know a priori which regime is relevant in a given case. One can 
simply compute v lc using both expressions given above and then take the larger value, which will 
always be the correct one. 

To determine whether direct observations of the IC emission is feasible, it is best to write u IC 
in terms of the underlying physical parameters describing GRB afterglow. The general expressions, 
though readily computable, are too complicated to give here; instead we show the results for two 
characteristic values of p. For p = 2.2 (and uniform-density ambient medium), we obtain 



TC 1 

v = max 



while for p = 2.4 we get 



TO 

v = max 



2.3xio^; 3 5 4>r 1 ^ 5 (i + -) a5 ; 

1.4xl0^ 5 - 2 1 ^ 3 5 VW- 3 (l + -)- 1 



QxlO^^r^a + z) - 4 ; 
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(5-2) 



Hz, 



(5-3) 



-1.3 



Hz. 



The first and second expressions in parentheses describe early and late time evolution of u IC , 
respectively. Note that when the latter expression applies, it may considerably overestimate the 
value of v IC , since it is based on an approximate description of the IC spectrum (see Fig. 1). 



-14- 



From Eqs. (5-2) and (5-3) it is clear that v reaches its minimum values at a time when 
v lc = . (Using Eq. (5-1) it can easily be shown that this is a general property of the transition 
frequency, as long as p > 2.2.) Since both Chandra and XMM cannot observe hard X-rays, the IC 
component can be seen directly only while, say, v lc < 5keV. This condition places a lower limit on 
the ambient density. For p = 2.2, we need n x > l.lE^ & e~^e^^{l + z) -0 ' 4 , while for P = 2A we 
should have n\ > 0.4£^°' 6 e~ ^e^ ^^ + z) -0-5 , in general agreement with Panaitescu & Kumar ( 
2000). 

If the ambient density is high enough to satisfy the above conditions then the following sequence 
of events should be observed in X-rays near 5keV. The X-ray flux begins to be dominated by IC 
emission at time t ~ 7.5 days e°; 9 5 e^ 2 n^ a7 (l + z) - 3 for p = 2.2 and t ~ 4.3 days e°; 9 5 e^ I 2 nj~ a7 (l + 
z) ' 3 for p = 2.4, causing the observed spectral slope to change from — p/2 to 1/3. After the 
transition, the IC flux increases linearly with time (Eq. 3-10), producing a bump in the X-ray 
lightcurve (see also figure 3 of Panaitescu & Kumar 2000) . 

The observed flux peaks when v 1 ^ moves into the X-ray band at time t ~ 7.2 days E^ef'^ x 

e B-2 n i ^ 9 (1 + z) 5 / 9 (independent of p). At the same time the spectral slope changes gradually 
from 1/3 to — (p — l)/2, which is still easily distinguishable from the characteristic —p/2 slope of 
the synchrotron high energy tail. From this point on, the observed flux evolves as 

~ 2 * 10 ~" (ssv) ~°' e ^^"i-^a + ' (5-4, 

for p = 2.2 and as 

f IC ~ 3 9 x irr 14 (— —\~°' 7 Rl-8,2.8 ,0.9 J.ln-2 (-, , ^1.6 f -1.3 er g (r v,\ 

for p = 2.4. For comparison, ACIS on Chandra can detect fluxes down to 4 x 10 _15 ergs _1 cm -2 
(in the energy band 0.4 — 6 keV) in a 10 4 s observation. Thus, for reasonable values of the afterglow 
parameters, the IC emission should stay visible for several weeks or even months, especially if E52 
is somewhat larger than unity, as is the case for a reasonable fraction of the bursts. 



6. Summary 

In this paper we calculate the detailed spectrum of the inverse Compton emission from a 
relativistic shock in the context of GRB afterglows. The general shape of this spectral component 
is very similar to the primary synchrotron spectrum. Like the latter, the spectrum of the IC 
emission can be roughly approximated by a broken power-law with break frequencies v]p , Vm , 
and u[ c (defined in §2.1). However, it differs from the synchrotron spectral component on three 
major points, summarized below. (I) At the low frequency end of the spectrum, v < u^ c ', the 
emission increases as oc is, rather than as oc v 2 characteristic of the synchrotron spectrum below 
v a . This may not have major practical consequences, however, since the low frequency part of the 
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IC spectral component is generally obscured by the synchrotron high energy tail. (II) The part 
of the IC spectrum which lies between v]P and vj^ extends over twice as many decades as the 
corresponding region (between v m and v c ) in the synchrotron spectrum. (Ill) The segments of the 
IC spectral component above the peak frequency (v^ in the slow cooling regime and above v]P 
in the fast cooling regime) cannot be well described by the pure power-law terms. The presence of 
additional logarithmic terms in this frequency range has the effect of considerably smoothing out 
the spectral breaks. Thus, no sharp changes in the spectral slope or in the rate of time evolution 
is expected during the observation of the IC emission above the peak frequency. In contrast, the 
observed transition between the synchrotron and Compton component can be sharp and produce 
quite dramatic spectral and temporal changes. 

In §3.1 we give a simple prescription for estimating the importance of IC emission for the 
total cooling rate of the shocked gas. We showed that if the fraction of energy contained in 
relativistic electrons, e e , is larger than that in magnetic field, e#, then during the fast cooling 
stage, IC emission is greater than the emission due to synchrotron by a factor of ^e e /eB- This 
factor decreases very slowly in the slow cooling regime, since at any given moment electrons located 
high enough in the power-law distribution are always cooling. Since the electron distribution is 
close to being flat, these electron can still radiate a significant fraction of the total electron energy. 
Thus, the afterglow evolution during the early part of the slow cooling regime is also dominated 
by IC emission. Moreover, we found that for e e /es > 10, the IC scattering is likely to remain the 
dominant cooling mechanism throughout the entire relativistic stage. 

As long as the IC emission dominates the total cooling rate, it sets the energy of the cooling 
break, j c , in the electron energy distribution, and determines the cooling frequency of the syn- 
chrotron spectral component. Therefore, the presence of IC emission changes the values of physical 
parameters inferred from current afterglow observations. In §4 we give the revised prescription for 
computing the total explosion energy, E, the fractions of energy in electrons and magnetic fields, 
e e and ee, and the ambient particle density, n, using the observed properties of the synchrotron 
spectral component. We show that at any instance in time two possible solutions for E, e e , e#, 
and n are allowed, one where IC cooling is unimportant and another where it dominates the total 
emission from the afterglow. There is no way to distinguish between the two solutions based on 
an instantaneous synchrotron spectrum alone. They differ only in their time evolution of emission 
above v c and in the strength of the IC spectral component. 

We also obtain an important constraint on the instantaneous values of the synchrotron break 
frequencies, peak flux of the synchrotron component and the redshift of the afterglow in the slow 
cooling regime. We show that due to the presence of the IC cooling rate, a combination of these 
parameters must always be smaller than 1/4. 

Finally, we discuss the possibility of detecting the IC emission component directly. We show 
that for reasonable values of the physical parameters, this component can be detected by Chandra 
a few days after the initial burst, as long as the ambient density is greater than ~ 1 cm -3 . Whether 
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or not the IC component is detected will be apparent from the change in observed spectral slope, 
as well as from the bump in the X-ray light curve. 

RS gratefully acknowledges support from the Sherman Fairchild foundations. AE was sup- 
ported by NASA through Chandra Postdoctoral Fellowship grant #PF8-10002 awarded by the 
Chandra X-Ray Center, which is operated by the SAO for NASA under contract NAS8-39073. 



A. Details of Compton Scattering 

The Thomson optical depth through the shocked medium, r ~ utuR, is generally very small, 
of order 10~ 6 . In this case the Compton y-parameter, y = ~f 2 T, determines the average fractional 
energy change of seed photons in each scattering by an electron with Lorenz factor j e . When y < 1, 
only single scattering needs to be considered in computing inverse Compton emission. When y > 1, 
multiple scattering can be important. However, in the context of GRB afterglows, a once-scattered 
synchrotron photon with initial energy hv generally has energy of order 7g hv > m e c 2 in the rest 
frame of the second scattering electron. Then Thomson limit no longer applies and the energy 
gain in each successive scattering will be reduced due to electron recoil and to the necessity of 
using Klein-Nishina scattering cross section. We conclude therefore, that multiple scattering of 
synchrotron photons can be ignored. 

For single scattering, the inverse Compton volume emissivity for a power-law distribution of 
scattering electrons is given by (Rybicki & Lightman 1979) 

poo rl 

£ C = 3a T d 7 N( 7 ) dxg(x)f Us (x), (Al) 

where x = v/A^ 2 v s , f Us is the incident specific flux at the shock front, and g{x) = l+x+2x In x— 2x 2 
takes care of the angular dependence of the scattering cross section in the limit 7> 1 (Blumenthal 
& Gould 1970). However, the quantities which are easiest to compare with the observations are 
the flux in the inverse Compton component, fl c = jl c ^R 3 /(4irD 2 ) and the synchrotron flux 
fu s = fv^R 2 /(4vrD 2 ), where R is the size of the shocked region and D is the distance to the 
observer. Strictly speaking, both j and / are measured in the rest frame of the shocked medium, 
but they transform in the same way to the observer's frame, so the transformation factors cancel 
each other out (as do redshift effects). Substituting these quantities into Eq. (Al) gives us 

fl c = Ra T d 7 N( 7 ) dxf Us (x), (A2) 

J 7 m JO 

where we approximated g{x) = 1 for < x < xq to simplify the integration. This simplified 
expression yields correct behavior for x < 1. The value of the parameter xo = V2/3 is set by 
ensuring energy conservation, i.e. setting x g(x)dx = Jq° xdx (we approximate xq ~ 0.5 in 
the main text of the paper). Note that the scattered flux distribution for monoenergetic photons 
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computed using this approximation has a maximum at x = Xo ps 0.47, very close to that of the 
exact distribution, which peaks at x « 0.61. 



A.l. Slow Cooling 

The distribution of seed photons is described by the synchrotron spectrum, which consists of 
four power-law segments (Sari et al. 1998). Then the inner integral in Eq. (A2) yields (to a leading 
order in v and zeroth order in v a jv m and v r n/v c ): 
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(A3) 



The quantity f max = fv s {v m ) is the flux value at the peak of the synchrotron spectral compo- 
nent. Note that the inverse Compton spectrum for monoenergetic electron scattering has the same 
frequency dependence above A^ 2 v a XQ as the input synchrotron spectrum above the self absorp- 
tion frequency, v a . In the range v < 47 2 f a xo, the functional form of the differential cross section 
dominates and I\ is linear in u, rather than quadratic. 

The integration over different electron energies again needs to be divided into four different 
regimes: 



fl c = Ra T x{ 



/; ( ,)^(7)V 



v < ui c - 
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(A4) 



where 7cr(i y ) = y/v/Avxo, and the break frequencies are defined as v[ c = Aj^VaXo, = ^Im^mXo, 
vl c = 47 2 i/ c x . 

Evaluating the integrals in Eq. (A4) and again keeping only the dominant terms, we obtain 
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There is no abrupt spectral slope change at v = y 7 v^u[ c , so that the region with the slope 
(1 — p)/2 extends over twice as many orders of magnitude as the corresponding region in the seed 
synchrotron spectrum. 

It is important to point out that for u > i/^p, the presence of logarithmic terms ensures that 
a broken power-law is no longer a good approximation to the correct spectrum. 

The value of the flux at the peak of the inverse Compton component is given by 

(p-l)(p+l/3) 



assuming that v a <C v m <C u c . 
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A.2. Fast Cooling 

In this regime, the inner integral in Eq. (A2) evaluates to: 
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Performing the appropriate integration over 7 we get the spectrum of the inverse Compton 
emission: 
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with the break frequencies defined as v a c = Aj^UaXo, v* c = 4j 2 u c xq, v£ = Aj^rnXo- 
In this regime, the flux at the peak of the inverse Compton component is 
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(A9) 



assuming that v a <C v c <C v n 



B. The Wind Case: p oc R~ 2 

For the instantaneous spectrum, there is no difference between the ambient medium with a 
wind-like or a constant density profile, since only the density in front of the shock at the time 
of observation matters. However, the two models clearly have different predictions for the time 
evolution of the observed emission. 

Eq. (3-1) and its exact (Eq. 3-2) and approximate (Eq. 3-3) solutions, as well as the expression 
for rj, are valid for any density profile, since they come from the analysis of a "snapshot" spectrum. 
However, since the time dependence of v c jv m is different for the wind-like and constant density 
profiles, the evolution of x is different in the two cases and so are the transition times between 
the three stages. For the wind density profile, the importance of IC emission decreases somewhat 
faster, since v c jv m oc t 2 in the wind case rather oc t as in the constant density case. This makes 
the difference between the intermediate regime (slow cooling, IC dominates) and the final regime 
(slow cooling with negligible IC) more significant. 

For the wind profile, the temporal evolution of x is given by: 

Assuming e e S> €b, slow cooling now lasts until 

ti c = 3.5(1 + z) (^) 3/2 (JjL) 1/2 A* days (B2) 
which is typically longer than in the constant density case. The IC dominated stage lasts until 
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This estimate is correct only during the relativistic stage and the importance of IC cooling declines 
fast, once the shock speed falls significantly below c. 

When inferring the parameters from a snapshot spectrum, the estimate of the total energy 
and the electron and magnetic energy fraction is not affected by changes in the density profile. 
Again one has to estimate the combination C as given by Eq. (4-6) and, provided that C < 1/4 
(otherwise there is no consistent solution), solve for x\ and X2, and substitute the values into Eq. 
(4-13). The equation for the ambient density also holds, but now it gives the density in front of the 
shock at the time of observation. However, since this density is not constant in time, the resulting 
number is not very useful. It is better to compute the normalization coefficient for the n oc R~ 2 
law. Using the Chevalier & Lee ( 1999) notation, p = 5 x 10 n gr cm~ 3 A*i?,^, we get 

10 " (dfcj'kpfcftwto* + *>< 1+ *>- <B4) 

which is similar to the expression of Chevalier & Lee ( 1999) with an extra factor (1 + x). 
For the new IC-dominated solution, given by X2 — 1/C 3> 1, we obtained: 

*= MS "- 1 (c^)" 1 (jspa)" 1 (ioSs)" (^)^ + -r 3 °i, 8 m 

As in the constant density case, this new solution, corresponds to a higher ambient density than 
the low-x solution. 
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Fig. 1. — Total energy spectrum of a GRB afterglow, calculated using the following parameters: 
p = 2.4, e e = 0.5, €b = 0.01, -E52 = 0.5, z = 0.5, n = 3. The synchrotron component is shown as 
a thin solid line and the inverse Compton component as a heavy solid line. A broken power-law 
approximation to the inverse Compton spectrum, normalized using Eq. (2-3), is plotted as a dashed 
line for comparison. Panel (a) shows the spectrum at t = 12 days, when the afterglow is in the slow 
cooling regime. Panel (b) shows the spectrum in the fast cooling regime, computed for t = 43 min. 




